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Even if dark matter particles are unambiguously discovered in experiments, there is no clear 
reason to expect that the dark matter problem has been solved. It is very easy to provide examples 
of dark matter scenarios (e.g. in supersymmetric models) where nearly identical detector signals 
correspond to extremely different relic densities. Therefore, the density of the particles discovered 
must be determined before their cosmological relevance is established. In this talk, I will present 
a general method to estimate the local density of dark matter particles using both dark matter 
and hadron collider experimental data when it becomes available. These results were obtained in 
collaboration with Gordon Kane at the University of Michigan. 



I. INTRODUCTION 

We are now confident that our universe contains 
a large amount of cold dark matter. The most pop- 
ular particle candidates for dark matter are weakly 
interacting massive particles (wimps). These parti- 
cles are being searched for directly and indirectly by 
dozens of experimental groups throughout the world. 
If we are fortunate, wimps may soon be discovered 
experimentally. 

Although the discovery of wimps in the galactic 
halo would have enormous implications for our un- 
derstanding of elementary particle physics, it implies 
very little about our understanding of dark matter. 
It it not reasonable to assume that the particle dis- 
covered represents all of the dark matter. This issue 
was first raised technically, though not resolved, in 
0| . Indeed, it is easy to give examples where par- 
ticles discovered in future and current experiments 
consist of less than 1% of the total dark matter. 

Even if weakly interacting massive particles are 
produced at colliders, it is still very important to 
directly determine their local density in the galactic 
halo. This can only be done with experimental data. 
During this talk, I will describe how this can be de- 
termined in the context of supersymmetry where the 
dark matter is the lightest supersymmetric particle 
(LSP). Similar analyses could be done for any dark 
matter candidate. 

This talk is organized as follows. First, I will il- 
lustrate why a discovery of dark matter particles in 
the halo is insufficient to address the dark matter 
problem and describe some of the uncertainties of 
relating the local and relic densities of dark matter. 
I will then describe how dark matter is detected di- 
rectly in experiments and offer the general form of 
the interaction rate. This will show what is required 
to determine the local density of wimps. I will give 
a very useful way to improve these expressions using 
data from different detector materials and at differ- 
ent recoil energies. 

In order to deduce the local density of wimps, it 
is essential to know the wimp's mass. I will present 
two methods to determine the mass of a wimp using 
direct detection data alone. The first reviews a well 
known relationship and the second presents prelim- 
inary results of the author. All of this work is then 



combined in the framework of the most general min- 
imally supersymmetric standard model where the 
neutralino is the wimp seen in dark matter experi- 
ments. A general procedure is presented to estimate 
the local density and explicit bounds are given. 

II. DISCOVERING (SOME OF?) THE 
DARK MATTER 

Let us imagine that a weakly interacting massive 
particle x nas been unambiguously observed in di- 
rect detection experiments. Such a discovery would 
be an enormous triumph of theoretical and exper- 
imental particle cosmology, have deep implications 
for our understanding of the universe, may herald 
the existence of supersymmetry, and will account for 
(at least) some of the dark matter in the universe. 
However, a discovery of dark matter particles is far 
from a solution to the dark matter problem: there 
is no reason to suspect that \ is a ^ the dark matter. 

What fraction of dark matter is represented by \ 
is a question that cannot be answered by experiment 
alone or theory alone. Furthermore, the answer will 
crucially depend on dark matter detection experi- 
ments. It is the purpose of this talk to describe how 
this question may be answered. 

It should be possible to determine the local den- 
sity of x using direct detection experiments. This 
is because, in a rough sense, they measure the lo- 
cal wimp density times its scattering cross section. 
Unfortunately, there are very few constraints on the 
scattering cross sections of most wimp candidates. 

However, it is not true that the signal rate depends 
on the cross section and density independently be- 
cause these are somewhat related. This is because 
the relic density f2 x is related to thermal produc- 
tion and freeze-out in the early universe. The rate 
of wimp-annihilation affects the relic density and de- 
pends on the wimp annihilation cross section, which 
in is turn somewhat related to the scattering cross 
section by crossing. 

Therefore, there is less freedom in the observed 
signal rate than one may have naively suspected. 
This can be illustrated qualitatively as follows. If the 
cross section is large, then most wimps would have 
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Direct Detection Signals and 
Relic Densities of 6050 cMSSMs 
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Direct Detection Signal in Ge at 20 keV (cpd/kg-keV) 

FIG. 1: The relic densities of 6050 constrained MSSMs 
as a function of direct detection signal strength in ger- 
manium. Experiments currently in planning or under 
construction may be able to observe signals of the order 
10" 4 cpd/kg-keV. 



annihilated in the early universe and local density 
would be small. Alternatively, if the cross section is 
small, then thermal freeze-out would have occurred 
very early and the density would be higher. In cither 
case, the cross section and density tend to compen- 
sate each other. 

The crude arguments above suggest that even a 
very small component of dark matter may be de- 
tectable because it may have a higher cross section. 
This has been referred to as the 'no-lose theorem' 
in recent conferences. Indeed, experimentalists may 
not lose out on discovering even a tiny fraction of 
the dark matter . 

This is seen in many realistic dark matter scenar- 
ios. In figure Q we have plotted the relic density 
against direct detection signal for some six thou- 
sand randomly generated, constrained minimally su- 
persymmetric standard models (without assuming 
any specific supersymmetry breaking scenario). By 
constrained, we mean that all of the models are al- 
lowed under current experimental constraints on su- 
persymmetry. These models were generated and an- 
alyzed using the DarkSUSY code 

Notice that for any particular signal strength, 
the relic density fluctuates over at least two orders 
of magnitude. In accordance with the 'no-lose the- 
orem,' experiments in the near future may detect 
even 1% of the dark matter or less. 

However, the no-lose theorem unfortunately also 
implies that the discovery of wimps in the galactic 
halo tells us very little about how much of the dark 
matter they compose. A wimp discovery could eas- 
ily represent a negligible fraction of the dark matter. 

Therefore, although wimps may be discovered in 
the near future, the dark matter problem will not 
be addressed until the density of wimps has been 



directly determined. 



III. LOCAL AND RELIC DENSITIES 

From studies of the cosmic microwave back- 
ground, large scale structure formation, and big 
bang nucleosynthesis, we know the cosmic-scale relic 
density of cold dark matter to be approximately 
^cdm^ 2 ~ 0.11 0. From our knowledge of the 
rotation of the Milky Way galaxy, the local dark 
matter halo density is known to be approximately 
Pcdm ~ 0.3 GeV/cm 3 .5j. It is obvious that any 
dark matter experiment on earth is only sensitive to 
the local density and not the relic density. 

Unfortunately, the relationship between the local 
and relic densities involves many details of galaxy 
formation and structure that are still not under- 
stood. Even if we were able to demonstrate that x 
has a local density of precisely 0.3 GeV/cm 3 , there 
remain important subtleties about our understand- 
ing of the dark matter in the universe as a whole. 

Noting these subtleties, however, it is extremely 
important to determine the local density of any 
wimp discovered in direct detection experiments. It 
would be very promising if the entire local dark mat- 
ter halo could be accounted for by wimps discovered 
in these experiments. 

Because direct detection experiments are sensitive 
to small-scale structure in the local halo density, a 
knowledge of the ambient halo density, p c dm ~ 0.3 
GeV/cm 3 may not be sufficient. Our knowledge of 
the local halo density is based on large scale sur- 
veys of star velocities in the Milky Way, and these 
measurements are not very sensitive to small-scale 
structure in the halo. 

There are several types of small-scale halo struc- 
ture which may effect direct detection experiments. 
For example, the earth may be within a stream of 
dark matter. This situation has been suggested by 
studies of the Sgr A stream and it has been estimated 
that our local halo density could be 0.3 — 23% higher 
than the ambient density |6|]. Alternatively, some 
authors have proposed that the halo may be clumpy 
or contain caustic structures 0. These small-scale 
perturbations in the dark matter density could have 
significant effects on direct detection rates. 

Most of the small-scale structure considered in the 
literature involves local, high-density regions of dark 
matter within the halo. Although these structures 
may make it easier to discover dark matter, they 
make it nearly impossible to assess what fraction of 
the halo is composed of x- 

Fortunately, there exist ways to check the smooth- 
ness of the halo profile. For example, dark matter 
streams or caustics may be identified or excluded us- 
ing directional dark matter experiments like DRIFT 
[8j- Also, a clumpy dark matter halo may possibly 
be identified by studying the time-dependence of a 
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wimp signal. Therefore, these questions may find 
answers in the foreseeable future. 

These ambiguities will need to be addressed be- 
fore the dark matter problem has been put to rest. 
However, for the purposes of this talk, we will as- 
sume that the halo is locally smooth and that if 
Px = Pedro = 0-3 GeV/cm 3 , then x is an °f the 
dark matter. 



IV. DARK MATTER DIRECT DETECTION 

Dark matter particles in the halo can be observed 
directly through their interaction with ordinary mat- 
ter hj. Although wimps interact only weakly, they 
will occasionally scatter off matter in detectors, de- 
positing a small amount of energy. Because wimps 
typically have masses on the order of a hundred GeV 
and move relatively slowly in the halo, on the or- 
der of a few hundred km/s, they typically deposit 
recoil energies from ~ 1 — 200 keV. Using very sen- 



sitive detectors, experiments can observe signals as 
low as a few keV. Using sophisticated coincidence al- 
gorithms, most experiments can remove virtually all 
background noise except scattering from neutrons. 

Direct detection experiments essentially measure 
the x-nucleus scattering rate as a function of recoil 
energy and time. In general, the signal rate is a func- 
tion of the cross section for %-nucleon scattering, the 
nuclear physics describing the nuclei in a particular 
detector, and the local velocity profile of the wimp 
fraction of the dark matter halo. 



A. Elastic Scattering Rate 

It will be helpful for us to state the explicit form 
of the differential interaction rate for a particular 
detector at the recoil energy q. Let the detector in 
question be composed of nuclei labeled by the index 
j, each with mass fraction Cj. Then, the differential 
rate of wimp scattering at recoil energy q is, 
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dv\F](q)[Z 3 f p + (A 3 -Z 3 )f n ] 2 + 



dQ Q=v 7rm x ^ y„ mIn . {q) v \JW* kj (2^ + 1) 

(1) 



where v m i nj (q) is the minimum velocity kinemati- 
cally capable of depositing energy q into the j th nu- 
cleus, f(v, t) is the local velocity distribution func- 
tion for the galactic halo, F?(q) and Sj mTl (q) are nu- 
clear form factors for coherent and incoherent scat- 
tering, respectively, Zj and Aj are atomic and mass 
numbers, Jj is the nuclear spin, a\ = a p + a n and 
ao = a p — a n , and the constant parameters f p<n 
and a Pt n describe the coherent and incoherent wimp- 
nucleon scattering cross sections, respectively |l4j . 
For a more detailed discussion of equation ^ please 
refer to any modern review of dark matter (see, e.g., 
i). 

It is important to note that equation ^ depends 
on several unknown parameters: 

1. the wimp's mass, m x , 

2. the particle physics of \ which determines the 
interaction parameters f Pt n and a Pt n, 

3. the velocity distribution of the halo, f(v, t), 

4. the local density of wimps, p x . 

These details will not be known when wimps are first 
discovered and may take many years to resolve. 



B. Prerequisites to Determine p x 

From the discussion above, it is clear that to de- 
termine the density p x , one must first 

1. Identify the particle x'i 

2. Determine m x ; 

3. Estimate the halo profile; 

4. Calculate the interaction parameters from the 
theory describing x- 

Each of these will require enormous efforts from both 
dark matter and collider physics experiments. 

The most important and perhaps most difficult re- 
quirement is the identification of x- This is not pos- 
sible from dark matter experiments alone. This is 
because these experiments observe only a few of x's 
quantum numbers. For example, it is unlikely that 
any amount of direct detection data can be used 
to differentiate between the lightest supersymmet- 
ric particle and the lightest Kaluza-Klein particle; 
if possible at all, this would probably require very 
precise data from several different nuclei. Therefore, 
although direct detection experiments may unequiv- 
ocally discover dark matter wimps, they cannot ex- 
plain the dark matter alone. 

Perhaps the most important parameter describing 
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X is its mass. This determines all of its kinematics 
and is crucially linked to the local density. Fortu- 
nately, m x may be calculable from direct detection 
experiments alone. The known methods of calculat- 
ing m x from dark matter experiments are described 
below. 

Because the mass may be observable, it may prove 
the key to the identification of %. If a neutral, stable, 
particle is observed at colliders with the same mass 
as that observed in dark matter experiments, then 
we could suspect that they are the same particle. 
Although this association is imprecise, it appears to 
be one of the best methods of identification. 

We should note, however, that determining the 
mass of x may be very difficult at hadron colliders. 
For example, if x is the LSP, it could take several 
years and an enormous effort to determine m x in a 
model-independent way. Most of the known tech- 
niques for determining the mass of the LSP rely on 
the framework of mSUGRA or a specific assump- 
tions about the relative masses of squarks and slep- 
tons. Therefore, it should be stressed that x may 
not be identified until long after it is discovered. 

Also, the halo profile must be known sufficiently 
well. As described above, any small-scale structure 
in the halo will dramatically alter the analysis of 
the local density. It is imperative that these issues 
be sufficiently resolved. 



Lastly, to compute the local density p x , one must 
know the interaction parameters f p . n and a p<n . To 
compute these parameters, one must know a great 
deal about the theory which describes x- It is clear 
that these parameters cannot be obtained from di- 
rect detection data alone: they rely on many pa- 
rameters of whatever extended standard model x is 
a part of. For example, if x is the LSP, then these 
parameters will be functions of the squark masses, 
mixing angles, gauge-content of the LSP, and higgs 
parameters. It is extremely unlikely that all of these 
will be known when x is discovered in direct detec- 
tion experiments. 



C. Combining Data 

All of the required analysis can be strengthened 
and empowered by combining data from different 
detectors over a range of recoil energies. There are 
many important insights and results based on the 
following framework. 

In general, the expression for the scattering rate, 
equation ^ is a second order polynomial in the four 
unknown interaction parameters f p , n and a p ^ n . To 
highlight this, it can be recast in the suggestive form, 
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It is clear from the expressions above, that by using 
data from 

1. different detector materials (varying the mass 
fractions, nuclear form factors, nuclear spins, 
and minimum velocities), 

2. different recoil energies (varying the nuclear 
form factors and minimum velocities), 

one can invert equation [5] to solve for ^fp^f p>n and 
y/p~^a p , n if the halo velocity distribution and m x 
were known. That is, given a halo model and wimp 
mass, one use find sufficient data from different de- 
tector materials and different recoil energies to de- 



termine ^/Px~f P ,n and y/p^cip^n (up to quadratic am- 
biguities). 

We should mention that there are many impor- 
tant situations in which the above analysis can be 
simplified. For example, because a PiTL are already 
scaled by linearly independent combinations of the 
incoherent nuclear form factors, Sj mn (q), it is not 
necessary to use different detector materials to solve 
for y/p^ap n. However, this will only work if there 
is data available from a detector with nuclei having 
non-zero spin which is sufficiently sensitive to inco- 
herent scattering. 

Although knowing the scaled interaction parame- 
ters will not directly determine the local density, it 
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can give enormous insight into the particle physics of 
X- For example, if x is the LSP, then the ratios a p / a n 
or dp/fn could possibly lead to important insights on 
tan /3, the degeneracy of the squark masses, mixing, 
and perhaps contain other information as well. This 
could be very important for collider physics and dis- 
entangling the MSSM. 



V. DETERMINING THE WIMP MASS 

Of all the factors required to interpret the ob- 
served signal rate, perhaps the most important is 
m x . Not only is the mass required to compute the 
density, but it also plays a critical role in the identifi- 
cation of x as described above. Fortunately, it may 
be possible to determine m x from direct detection 
experiments alone. 

There are roughly two ways to determine m x from 
direct detection data. One method, that using the 
annual modulation crossing energy, was described 
in the dark matter review article by Primack et. 
al. in 1988 @. Although it seems unlikely to have 
originated in a review article, we have been unable 
to find any previous author mentioning this effect. 
The other method has been developed by the author 
in collaboration with Gordon Kane and represents 
work still in preparation. Both of these methods 
rely only on the kinematics of the halo, although 
neither are particularly sensitive to the precise halo 
model ^3 (although, see caveats in [Io|). 



The Difference Between Summer and Winter 
Differential Wimp Interaction Rates 
For Various Energy Bins 
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FIG. 2: The difference between direct detection signals 
in germanium in June and December as a function of 
recoil energy. This plot was generated for an MSSM 
with m x ~ 161 GeV. 



small-scale halo structure. Furthermore, if the wimp 
is very light, then the the crossing energy may be 
well below the threshold of the detector and there- 
fore not observed at all. 



A. Annual Modulation Crossing Energy 

As the earth orbits the sun, its velocity through 
the galactic dark matter halo varies between roughly 
250 and 190 km/s This in turn causes annual 
modulation in the scattering rate. However, the am- 
plitude of this modulation varies as a function of re- 
coil energy and changes sign. For a more detailed 
description, see, for example, |10|. 

In figure |2 we plot the difference between scat- 
tering rate in June and December as a function of 
recoil energy for several detector materials. Notice 
that there is a particular energy, called the 'crossing 
energy,' at which no annual modulation is observed. 
As pointed out by Primack et. al. 0, the crossing 
energy is an explicit function of the masses of the 
wimp and detector nuclei which can be derived easily 
from kinematics. Therefore, if crossing is observed, 
one can determine the mass of the wimp explicitly. 

This method is moderately robust. Specifically, if 
there is an energy at which the annual modulation 
amplitude changes sign, then one can confidently 
determine the wimp mass to within approximately 
10%. However there are some important subtleties 
and caveats to this analysis described in, for exam- 
ple, [13. These include the effects of bin-sizes and 



B. Kinematical Consistency 

Recall that if the halo velocity profile and m x are 
known, then direct detection data from different de- 
tector materials and different energies can be used 
to solve for y/p^f P , n and ^p^a p ^ n . If the halo veloc- 
ity distribution is known, then only the wimp mass 
is required to determine these. 

Let us assume that the local halo velocity profile 
can be adequately approximated and that there ex- 
ists enough data to solve for ^fp^fp.n and y/p^a p . n if 
the mass were known. (If, for example, there is only 
data sufficient to solve for a p Tl , it will be clear how 
to proceed along similar lines). Because many direct 
detection experiments observe scattering rates in a 
large number of recoil energy bins, we can generally 
expect to have many more measurements than the 
minimum required to solve the system of equations. 

Because the interaction parameters are absolute 
constants, all minimal, linearly independent combi- 
nation of measurements used to solve for the scaled 
interaction parameters will agree if the correct mass 
were used in the derivation. However, if an an ar- 
bitrary m' x were used to solve for these parameters, 
different calculations will not in general agree. 

This motivates us to define a 'kinematical consis- 
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Algorithm to Determine the Mass of the LSP 
Given Direct Detection Kates 
for both Nal and Ge Detectors 

(model ATLAS_SDSY02) 




112 114 US 118 

Mass of the. LSP used for Calculation 



FIG. 3: The function C( m x) where the wimp corresponds 
to the neutralino in the MSSM specified by ATLAS 
SUSY point 2 11]. The models and data were gener- 
ated within the framework of the DarkSUSY package, 



tency' function, £(m' x ), which compares the values of 



\fPxfp,nr\fPxO,p tn obtained using different indepen- 
dent subsets of the data as a function of m' x used. 
Specifically, let C( m ' x ) be given by 



C(m' x ) = |(a p (i) — a p (j)) 2 + similar termsj 



where the indices i,j represent a minimal set of data 
used to compute the constants given the particular 
m' x . It is necessary that C{ m ' x ) — when m' = m x , 
but this is not a sufficient condition. Specifically, we 
have not found any way to demonstrate that m x is 
the unique root of £(m' x ), although we have found 
no example where it has multiple roots. 

To determine the wimp mass, one varies m' x un- 
til C( TO x) = 0. To test how useful this technique 
is, we applied it to some six thousand random, con- 
strained MSSMs. In every single model tested, the 
correct mass was determined to near-arbitrary pre- 
cision. Figure 01 illustrates a typical plot of ((m' x ). 
Notice that £ has an extremely sharp minimum, de- 
creasing many orders of magnitude within a few GeV 
of the true mass of the LSP. It should be stressed, 
however, that experimental uncertainties and reso- 
lutions were not considered for these calculations. 

Although this method appears quite promising, 
these results should be considered preliminary. Sev- 
eral questions still remain regarding how robust this 
calculation is after the introduction of experimental 



uncertainties and halo model ambiguities. Without 
these uncertainties, the algorithm yields the correct 
mass to seemingly arbitrary precision. It will be in- 
teresting to see how this changes in more realistic 
circumstances. 



VI. BOUNDS ON THE LOCAL DENSITY 

Because m x can be determined in principle using 
the methods described above and we are working 
under the assumption that the local halo velocity 
profile can be adequately approximated, we can de- 
termine y/px/p.n and/or ^p^a Ptn independent of the 
identification of x- Therefore, to determine p x it is 
sufficient to know any one of the interaction param- 
eters. This is an enormous improvement over the 
general case, for which all of the interaction param- 
eters were required. 

Therefore, any bounds on the interaction param- 
eters will translate into bounds on the local density. 
Unfortunately, these parameters can only be com- 
puted in the framework of a very explicit model for 
the wimp. Furthermore, these parameters are typi- 
cally very poorly constrained. 

In order to address the question of the local den- 
sity p x , one must specialize to a particular candi- 
date particle in a specific extension of the standard 
model. Therefore, we cannot proceed without losing 
some generality. 



VII. NEUTRALINO DARK MATTER 

The most popular and perhaps best-motivated 
candidate for cold dark matter is the lightest su- 
persymmetric particle (LSP), predicted supersym- 
metric extensions of the standard model which con- 
serve i?-parity. Indeed, supersymmetric dark matter 
was predicted before it was known that non-baryonic 
dark matter was needed. In most MSSMs allowed by 
experimental constraints, the LSP is the neutralino, 
X, which is the supersymmetric partner of the neu- 
tral gauge and higgs bosons. For an extensive and 
authoritative review of supersymmetric dark matter, 
see Jungman et. at. 



A. Interaction Parameters 

Given a completely specified supersymmetric 
standard model, one can straightforwardly compute 
the interaction parameters. It should be noted, 
however, that one does not directly compute / Pjn 
or a p . n . Rather, %-quark (an %-gluon) interaction 
parameters are calculated and these are used to 
determine the x-nucleon parameters. 

To tree level, the %-quark interaction parameters 
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include the following diagrams. 



x q_i q% x 




X Qi x ft 

X ft x ft 




x ft x ft 



It is clear that these will depend on many of the pa- 
rameters in the model. Specifically, they are func- 
tions of the 



1. gauge content of the lightest neutralino, 

2. most of the squark masses and mixing angles, 

3. tan/3, the ratio of the vacuum expectation 
value of the two higgs bosons, 

4. higgs mass parameters (only for the coherent 
interactions). 

It must be emphasized that most of these param- 
eters will be extraordinarily difficult to measure in 
practice (especially at hadron colliders). There do 
not exist today general, model-independent meth- 
ods of determining most of these parameters. 

Although we will not derive these here, to illus- 
trate the dependence on each of these parameters, 
the incoherent scattering of x with a u-quark is given 
by, 
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In the expression above, the matrices Hl,r are 
3x6 projection matrices given in the basis 
(ul, CL,tL, ur, cr, tfi); @ u is a unitary matrix which 
diagonalizes M 2 so that Af 2 dia s = OiM 2 @ u pj|; 
the subscript j on q~j corresponds to the flavor and 
handedness of the quarks so that j = 1, . . . , 6 corre- 
sponds to (ul, Cl, II, ur, cr, tii); and gauge content 
of x is given by 

X = N § \B) +N VV \W)+ N 6l \H!) + N S JH 2 ). 

A similar expression describes scattering with d, s- 
quarks. 

The coherent parameters are similar to the inco- 
herent ones except that they also contain higgs ex- 
change at tree-level. This implies that in addition to 
squark masses, mixing angles, tan/3, and the gauge 
content of %, one must also know the higgs masses. 
Therefore, in general, less knowledge of the MSSM 
is required to compute a p<n than f PtTl . 



B. Limits on Scattering Parameters 

As we have shown, the interaction parameters de- 
pend on very detailed knowledge of the MSSM. Un- 
fortunately, these may not be known until well after 
dark matter particles have been observed directly in 
experiments. We should somehow try to estimate 
them using partial information and any available 
bounds on the MSSM. 

It should be clear that even if all of the squark 
masses and mixing angles are unknown, we can still 
place limits on the interaction parameters using ex- 
clusion bounds. In general, one can typically find 
a way to make use of what is known and constrain 
what is not known to estimate and limit the inter- 
action parameters. 

Beginning without almost any parameters of the 
MSSM determined, we find that we can still place 
rather strong limits on a p .„. For example, we have 
found that given only upper and lower bounds on 
tan j3 and a lower bound on the lightest squark mass, 
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there is a strict upper bound for the incoherent \- 
quark scattering parameters. If there is a strict lower 
bound on the lightest squark mass, say and 
tan/3 is bounded so that sin/3 > sin/% [l7|. then 
there is a strict upper bound on a p , n . It should be 



mentioned that these types of bounds already exist 
today, at least in the framework of particular super- 
symmetry breaking scenarios. In this case, it can be 
shown, that the magnitude of a u is strictly bounded 
by 



a u < 



16m 



— (N ~ - N - ) + 

2 w Hl H2 8 (m? -(m2+m u )2 



^- tan 2 6 W N% + \n%, H — ^ „ 

18 B 2 w m^sin 2 ^ H * 



-N% 



[ 717, 

+ -tanO w \Ng\\Nw\ cos(a^) + - — - |V # | |A^ | cos(o^ 2 - a^) 
o m\y sm pi 



niw sin fa 



tan 6 W | Ng 1 1 Nfr \ cos(a^ ) >, 



where a ^ 2 , and are the relative phases between 
Njj , Ny^ and A 7 ^ , respectively. 

This expression has six real unknowns. Notice 
that by the normalization of the neutralino wave 
function, 



WB? + Wy r ? + \N Al \ i + \N^ = \, 



the parameter space is compact. Therefore, a u can 
be absolutely maximized with respect to all six un- 
knowns. Specifically, although all of the gauge con- 
tent of the neutralino may be unknown, one can ab- 
solutely limit the x - quark and hence the x-nucleon 
interaction parameters. 

It should be emphasized that the analysis used 
to derive the above bound was for the most gen- 
eral softly-broken supersymmetric standard model; 
no ad hoc supersymmetry breaking scenarios such 
as mSUGRA were assumed. It is obvious that if a 
particular supersymmetry breaking scenario were as- 
sumed, the above expressions would be enormously 
simplified. However, these types of assumptions are 
very difficult to justify (theoretically or experimen- 
tally) and therefore greatly limit the generality of 
the work. 

It is important to note the flexibility of the deriva- 
tion involved to compute these bounds. If, for exam- 
ple, the masses of several light squarks were known, 
one can greatly improve the above bounds by in- 
cluding these in the explicit expression for a q and 
then maximizing it relative to the parameters that 
remain unknown. In this manner almost any addi- 
tional knowledge can be added to arrive at stronger 
statements. Therefore, not only do these bounds 
grow more restrictive with increasing knowledge of 
the MSSM, but they continue to approach a realistic 
estimate of the interaction parameters. 



C. Strong Lower Bound on p x 

From the work above, it is clear that given ade- 
quate bounds on tan/3 and a lower bound on the 
lightest squark mass, there exist strong, model- 
independent upper bounds on a Pi „. These in turn 
can be used to place a very strong lower bound 
on the neutralino relic density because we know 
y/Px a p,n ■ 

To test this method, we considered some six thou- 
sand randomly generated MSSMs that are consistent 
with all known bounds on supersymmetry. For each 
of these models, we calculated the interaction rates 
for a Nal detector in twelve recoil energy bins. This 
(idealized) data was used to compute the mass of 
the LSP, using the kincmatical consistency function, 
and to solve for y/p^a PiTl . Upper bounds were calcu- 
lated for a Pi n assuming 10% uncertainty in tan fi and 
a lower bound on the lowest squark mass of either 
200 GeV or the actual mass of the lightest squark, 
whichever is less. The specific gauge content of the 
neutralino was taken to be known for each model, 
however, for the sake of computational simplicity 
[l8|. Using the upper bounds for a p ,„, we obtained 
a lower bound on the local density p x . 

Figure 0] illustrates the results of using this algo- 
rithm for each of the randomly generated MSSMs. 
Notice that the estimated local density is always 
strictly less than the true local density — as required 
by it being a lower bound. Also notice that for many 
models the lower bound was not such a poor esti- 
mate. This will be the case, for example, when the 
lightest squark mass is near or below the 200 GeV 
bound. 



VIII. CONCLUSIONS 

We have seen that, by itself, a discovery of dark 
matter particles in our galactic halo cannot address 
the dark matter problem of the universe. However, 
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Lower Bound and Estimate 
of the Local Neutralino Density Using Nal Data 
with Bounds on Lightest Squark Mass and tan/3 




True Local Density (GeKL-ai-^ 



FIG. 4: This plot compares the lower bound and es- 
timate of the local denisty computed using the strong 
upper bound for a P: „ to the true local density for each 
model. The red line indicates perfect agreement. Notice 
that the procedure correctly determined a lower bound 
for the local density for every model. 



combined with data from colliders to identify the 
particle and limit its interaction parameters, we can 
generally estimate its local density. 

In the framework of supersymmetry, we have 
shown a robust and iteratively improvable method 
to estimate the local density of a neutralino LSP 
observed in direct detection experiments using any 
information available about the MSSM. 

Therefore, although the dark matter problem 
may not be solved immediately when wimps are 
observed, there are clear and general ways to 
address their cosmological significance. 
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